Measuring infrared contributions to the QCD pressure* 
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For the pressure (or free energy) of QCD, four-dimensional (4d) lattice data is available at zero baryon density 
up to a few times the critical temperature T c . Perturbation theory, on the other hand, has serious convergence 
problems even at very high temperatures. In a combined analytical and three-dimensional (3d) lattice method, 
we show that it is possible to compute the QCD pressure from about 2T C to infinity. The numerical accuracy is 
good enough to resolve in principle, e.g., logarithmic contributions related to 4- loop perturbation theory. 



Introduction. The properties of QCD matter 
are expected to change above a critical tempera- 
ture T c ~ 200 MeV. While the low-temperature 
phase is governed by bound states, such as 
mesons, the high-temperature phase should, due 
to asymptotic freedom, look more like a gas of 
free quarks and gluons. Any observable witness- 
ing this change is therefore a potential candidate 
for direct or indirect measurements in heavy-ion 
collision experiments. One such observable is 
the free energy density, or pressure of the sys- 
tem, since according to the Stcfan-Boltzmann law 
(i.e. neglecting interactions), it is proportional to 
the number of effective degrees of freedom. 

For vanishing baryon density, at a temperature 
T and volume V, the free energy density / is sim- 
ply given by the functional integral 

/ = ~ In/ V[A^} e -l> qT fv d3xC ^ A ^ , (l) 

where Ce is the standard QCD Lagrangian. For 
V — » oo, the pressure p is given by p = — f. 

The most direct way to evaluate Eq. (|l|) is to 
measure it numerically on the lattice. This has 
been done by a number of groups, for zero as 
well as non-zero number of fermion flavours Nf 
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(e.g. HQ). The general picture emerging is the 
following: The pressure rises sharply in the in- 
terval (1 — 2)T C , to level off at a few times T c . 
At the highest temperatures used in the simu- 
lations, typically (4 — 5)T C , the deviation from 
the Stefan-Boltzmann limit is about 15%. At 
even higher temperatures, the pressure is then ex- 
pected to asymptotically approach the ideal-gas 
value p (T) = (tt 2 T 4 /45) (JV C 2 - 1 + (7/A)N c N f ), 
where N c denotes the number of colours. 

The 15% deviation cannot be systematically 
understood in terms of (finite-temperature) per- 
turbation theory While the expansion is known 
analytically to 5th order in the gauge coupling 
g ||, convergence properties are extremely poor. 
Therefore, in the past few years a lot of effort 
has gone into refined and/or alternative analytic 
approaches (5| . A general feature of these works 
is the suppression of infrared effects. While this 
suppression does not seem to be crucial in the 
computation of the pressure, one might not be 
satisfied with accepting it as an ad hoc assump- 
tion. The aim of this talk is to briefly review our 
framework of resumming the infrared contribu- 
tions to the pressure to all orders || . 

Pressure via effective theory. A way to 

understand the poor convergence of the ordi- 
nary perturbative expansion is the observation 
that when gT <C ttT, the system undergoes di- 
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mensional reduction (see, e.g., || and references 
therein). In the case of QCD, the effective theory 
is a 3d SU(-/V C ) + adjoint Higgs model: 

£ 3d =\F^ +\[D l , A f +\ m \Al +\\,A\. (2) 

The parameters (173, m§, A3) are related to the 
physical parameters of the full 4d theory (T, Aj^t) 
by perturbatively integrating out the hard modes 
(~ ttT). Below, we will work with the dimen- 
sionless variables x = {\%/g\) and y — (m^/g^), 
known at next-to-leading order JjJ . 

The effective theory is confining, hence non- 
perturbative ||. Therefore, the only way to sys- 
tematically include the infrared contributions to 
the pressure is to treat £3d on the lattice. 

Let us rewrite the pressure, up to hard-scale g 6 
contributions, as 
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where 5 ~ 10~ 4 , Po(T) is the Stefan-Boltzmann 
value, and the dependence on the scale /13d, which 
originates from an infrared divergence of the 4d 
part, cancels against an ultraviolet tcrrr^] in the 
dimensionless 3d free energy density J^g(x,y). 
We write J-jj^(x, y) as 
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and aim at measuring it on the lattice. This re- 
quires a determination of the quadratic and quar- 
tic Higgs field condensates (which equal the par- 
tial derivatives under the integral), as well as a 
(4-loop) perturbative computation in lattice regu- 
larization, to match to the MS scheme. We choose 
to fix the integration constant •^^jgC^Oj 2/o) at high 
temperatures T ~ 10 T c , where one is confident 
that continuum perturbation theory converges^]. 

As a first result obtained along the strategy 
outlined, Fig. [j] shows the normalized pressure. 

1 This is precisely the way the effective theory is set up: 
dependence on a matching scale has to cancel. 
2 T c is of the order of A^jjg-, with a coefficient measurable 
on 4d lattices. 
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Figure 1. The pressure after the inclusion of the 
infrared part according to Eq. (H) (from Ref. @). 
Statistical errors are shown only for eo = 10. 

The integration constant has been fixed pertur- 
batively on the 3-loop level, allowing for an ad- 
ditional constant eo, which represents an (up to 
now) unknown g| contribution. In principle, this 
constant can be determined in a setup equivalent 
to the above, after splitting off its perturbative 
part: A further reduction relates eo to the free 
energy of 3d pure gauge theory (see, e.g., ||), 
which can be determined on the lattice. On the 
lattice side, we have only included the quadratic 
scalar condensate in Fig. |l|, due to reasons ex- 
plained in the next section; however, at T — > T c 
the quartic one will become important as well. 

3d lattice measurements. Let us now discuss 
in more detail how the actual measurement of the 
3d condensates is carried out. We need to relate 
lattice and MS regularization schemes. Here, the 
super-renormalizability of the 3d theory plays an 
important role: there are only a finite number of 
divergent terms, such that analytic relations can 
be computed exactly near the continuum limit Q . 
Writing (3q — (2N c /g^a) where a is the lattice 
spacing, the observables we need behave schemat- 
ically like 



8 T — . 



0G 



lim {(TrA2) L +/3 G +ln/3 G +l}, (5) 



3 



y = 9.56 



1 - 



5 



-2 - 



-3i 



G--o bare / 100 

u— a p o ! / 10 



all'Z 
I • ( ) - perlurbative 



0.1 



0.2 0.3 

2 



0.4 



0.5 



Figure 2. A measurement of one of the partial 
derivatives needed in Eq. (||), d x F, at a fixed y, 
as a function of the lattice spacing g\a. Note that 
the raw data (open circles) is scaled by a factor 
1/100. The other curves result from subtracting 
various divergences, Eq. (Q); taking into account 
improvement |l^|; and finally (filled circles), sub- 
tracting the perturbative value. 
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For Eq. (||) all three coefficients are known JlO[ . 
In the case of the quartic condensate, Eq. (g), all 
we know so far are the four divergent terms: the 
weakest logarithmic divergence can be obtained 
in a 4-loop continuum calculation, but the 4-loop 
lattice constant is still unknown. This is precisely 
the reason why in Fig. [IJonly the effect of (TrAp) L 
is included. While the graphs needed can be sys- 
tematically generated |ll|, carrying out the lat- 
tice integrals remains a major challenge. 

Nevertheless, we can show how the data involv- 
ing the quartic condensate behaves. In Fig. |^, we 
demonstrate for one specific temperature that the 
series in Eq. (g) indeed regulates the lattice op- 
erator as g\a — > 0. Indeed, taking into account 
successive terms of the series, one observes how 
the divergence at small a becomes weaker. Note 
that while the continuum limit can be taken al- 



ready now, the actual value of d x will have to 
await the above-mentioned 4-loop computation of 
the finite term in lattice regularization. 

Conclusions. We wish to point out two trends 
seen in Fig. |l| First, the outcome is sensitive 
to the non-perturbative parameter eo, which in 
principle can be determined by additional com- 
putations. Clearly, there exists a range for that 
parameter which leads to a sensible result. 

Second, comparing with Fig. 1 of g, at T > 
30A^g the curves for g 5 (denoted by there) 
and eo = 10 fall almost on top of each other, 
signalling a cancellation of all higher-order terms 
(determined by the quadratic Higgs condensate) 
against the large non-perturbative g e contribu- 
tion. Hence, in this temperature range the pres- 
sure is indeed dominated by ultraviolet effects. 

Let us also remark that as demonstrated in 
Fig. 3d lattice results involve always a system- 
atic extrapolation to the continuum limit, and are 
precise enough to resolve, e.g., logarithmic effects 
related to 4-loop perturbation theory. 

We end by noting that the inclusion of Nf 
fermion flavours as well as a baryon chemical 
potential /ib pose no further complications, and 
hence provide for a natural extension of this in- 
vestigation. 
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